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LONG-TERM  GOALS 

To  understand  the  dynamics  of  physical  oceanographic  circulation  processes  on  continental  shelves  and 
slopes  with  emphasis  on  the  mechanisms  involved  in  across-shelf  transport. 

OBJECTIVES 

To  apply  numerical  circulation  models  to  process  studies  and  to  simulations  of  continental  shelf  and 
slope  flow  fields,  including  the  inner  shelf  region  and  the  nearshore  surf  zone,  to  help  achieve 
understanding  of  the  flow  dynamics. 

APPROACH 

Numerical  finite-difierence  models  based  on  the  primitive  equations  and  the  shallow-water  equations  are 
applied  to  two-  and  three-dimensional  flow  problems  relevant  to  the  dynamics  of  continental  shelf  and 
slope  flow  fields  and  to  the  circulation  in  partially  enclosed  seas.  At  present,  the  Princeton  Ocean 
Model  (POM)  (Blumberg  and  Mellor,  1987)  is  being  utilized  for  studies  with  the  primitive  equations.  A 
shallow-water  equation  model  has  been  developed  and  applied  to  studies  involving  vorticity  dynamics  of 
wave-averaged  currents  in  the  nearshore  surf  zone.  The  numerical  experiments  are  supplemented  with 
analytical  studies  whenever  possible. 

WORK  COMPLETED 

Numerical  experiments  addressing  the  nature  of  nonlinear,  finite  amplitude  shear  instabilities  of 
alongshore  currents  in  the  surf  zone  through  the  solution  to  idealized,  forced,  dissipative,  initial-value 
problems  have  been  performed.  Bottom  topographies  with  constant  slope,  with  an  alongshore-uniform 
sand  bar,  and  with  an  alongshore  variable  sand  bar,  have  been  used  with  periodic  boundary  conditions  in 
the  alongshore  direction.  Forcing  effects  from  breaking  surface  waves  are  calculated  from  gradients  in 
the  radiation  stress  tensor.  These  are  determined  using  the  formulation  of  Thornton  and  Guza  (1986) 
generalized,  for  the  alongshore  variable  topography,  to  two-dimensional  situations.  Dissipative  effects 
are  modeled  by  linear  bottom  friction. 

In  coordination  with  closely  related  work  in  the  NOPP  project  "Development  and  Verification  of  a 
Comprehensive  Community  Model  for  Physical  Processes  in  the  Nearshore  Ocean",  the  Princeton 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

SEP  2000  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2000  to  00-00-2000 

4.  TITLE  AND  SUBTITLE 

Modeling  of  Coastal  Ocean  Flow  Fields 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

College  of  Oceanic  and  Atmospheric  Sciences„Oregon  State 

University, 104  Ocean  Admin  Building„Corvallis„OR, 97331 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

7 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


Ocean  Model  (POM)  has  been  adapted  for  studies  of  the  three-dimensional,  short-wave-averaged 
cireulation  in  the  surf  zone.  Parameterized  forcing  from  breaking  waves,  represented  by  gradients  in  the 
radiation  stress  tensor,  and  parameterized  effects  of  near-surface  wave-indueed  mass  flux  are 
incorporated.  Long  range  objeetives  include  development  of  a  single,  unified  modeling  capability  for 
the  low-frequency  circulation  in  the  surf  zone,  foreed  by  breaking  surface  gravity  waves,  and  for  the 
cireulation  over  the  adjaeent  inner  shelf,  foreed  by  wind  stress. 

Researeh  has  been  initiated  on  modeling  studies  for  inner  shelf  flow  fields  by  considering  the  wind- 
forced  circulation  off  Duck,  NC  in  August-November  1994  during  the  time  of  the  Coastal  Oeean 
Processes  (CoOP)  field  experiment.  This  is  part  of  Ph.D.  thesis  research  by  Brandy  T.  Kuebel.  For  the 
initial  numerieal  experiments,  the  assumption  of  alongshore  uniform  two-dimensional  flows,  with  spatial 
variations  in  the  aeross-shelf  (x)  and  vertical  (z)  directions,  is  utilized.  The  model  is  forced  by  observed 
wind  stress  and  heat  flux.  Modekdata  comparisons  are  made  for  velocity  and  temperature  fields. 

During  the  time  period  of  the  field  experiment  and  the  model  calculations,  both  stratified  (August)  and 
unstratified  (October-November)  conditions  exist  allowing  eomparison  of  the  shelf  flow  response  in 
these  two  different  regimes. 

In  a  separate  effort,  numerical  model  experiments  utilizing  POM  have  been  conducted  to  study  the 
mesoscale  circulation  in  the  Gulf  of  California.  This  is  Ph.D.  thesis  research  by  Antonio  Martinez.  The 
separate  effects  of  forcing  by  winds,  by  coastal-trapped  waves  ineident  from  the  south,  and  by  the  tides 
have  been  examined.  A  relatively  high  resolution  grid  (3  km  horizontal  grid  size,  50  sigma  levels  in  the 
vertical)  has  been  employed  to  adequately  resolve  the  mesoscale  flow.  The  wind  forcing  experiments 
have  been  run  for  365  days  (August  1996  -  July  1997)  after  a  60  day  spin-up  period.  The  wind  stress  is 
obtained  from  a  combined  produet  of  scattero meter  measurements  and  NCEP  analyses  (Milliff  et  al, 
1999).  The  eoastal-trapped  wave  experiments  have  been  run  for  an  80  day  period  1  July  -  19  September 
1984  during  whieh  time  extensive  current  (Merrifield  and  Winant,  1989)  and  hydrographie  (Bray, 

1988)  measurements  were  made  in  the  gulf  These  measurements  are  utilized  for  modekdata 
comparisons.  It  has  been  concluded  from  previous  observations  (Merrifield  and  Winant,  1989)  that 
storm-generated  ineident  coastal-trapped  waves  make  a  major  contribution  to  mesoscale  variability  in 
the  gulf  The  ineident  coastal-trapped  waves  are  assumed  to  have  the  spatial  structure  of  the  first 
baroclinie  linear  mode  with  time  variability  given  by  coastal  sea  level  measurements  at  Aeapuleo.  The 
emphasis  in  the  experiments  involving  tidal  forcing  has  been  on  a  study  of  the  eharacteristies  of  the 
internal  tides  and  of  the  assoeiated  tidal  mixing  in  the  gulf 

RESULTS 

Results  from  the  study  of  nonlinear  shear  instabilities  of  alongshore  currents  in  the  surf  zone  over  beaeh 
topography  with  alongshore-variable  sand  bars  (Slinn,  Allen,  and  Holman,  2000)  show  significant 
influence  of  alongshore  topographie  variability  on  the  nearshore  eirculation.  In  particular,  one  notable 
feature  is  the  tendeney  for  contours  of  both  the  time  mean  and  the  root  mean  square  vorticity  fields  to 
align  along  eontours  of  constant  depth.  This  tendency  is  attributed  to  both  the  direet  effect  of  the 
variable  topography  on  the  wave  forcing  and  also  to  the  dynamieal  response  of  the  currents  to  variable 
topography.  The  contributions  from  each  of  these  two  mechanisms  are  of  eomparable  magnitudes. 


Initial  applications  of  POM  to  surf  zone  circulation  studies  have  focused  on  alongsho re-uniform  two- 
dimensional  motion.  The  model  has  been  applied  to  studies  of  the  eirculation  off  Duck,  NC  with  model 


results  being  eompared  to  veloeity  measurements  from  the  DUCK94  field  experiment  (Gareez  Faria  et 
al,  1998,  2000).  The  POM  results  determine  model  solutions  for  the  (x,z)  structure  of  the  wave- 
averaged  alongshore  and  across-shore  velocity  fields.  Reasonable  agreement  is  found  between  the 
model  and  the  measured  velocities.  The  effects  of  tidal  elevation  change  on  the  circulation  are 
investigated  and  show,  in  particular,  variations  in  the  strength  of  the  undertow  over  the  bar  and  in  the 
trough  with  tidal  height  that  are  in  general  agreement  with  velocity  measurements  from  the  fixed  array. 

Numerical  model  studies  of  the  two-dimensional  circulation  off  Duck,  NC  during  August-November 
1994  provide  detailed  information  about  wind- forced  inner  shelf  flow  fields.  A  comparison  between 
measured  and  modeled  alongshore  currents  at  two  across-shelf  positions  shows  a  good  correlation  (see 
Figure  1).  Also  plotted  in  Figure  1  are  monthly  means  and  standard  deviations  of  the  model  potential 
density,  streamfimction  and  alongshore  velocity  fields,  which  are  significantly  different  during  stratified 
and  unstratified  conditions.  In  August  (stratified),  both  upwelling-  and  downwelling-favorable  winds 
force  the  stratified  flow  on  the  shelf,  with  an  apparent  mean  upwelling  circulation.  A  well- formed 
northward  coastal  jet  develops,  with  opposing  southward  velocities  both  below  and  offshore  of  the  jet. 
In  October  (unstratified),  strong  downwelling-favorable  winds  prevail  and  the  shelf  becomes  almost 
completely  mixed  in  the  vertical.  Alongshelf  velocities  are  southward  everywhere  and  decrease  in 
magnitude  with  depth. 

Results  from  the  modeling  studies  of  the  mesoscale  circulation  in  the  Gulf  of  California,  categorized  by 
forcing  mechanism,  are  briefly  summarized  below. 

Coastal-Trapped  Waves:  These  waves  propagate  northward  into  the  gulf  along  the  coast  of  mainland 
Mexico.  Typical  periods  of  the  energetic  events  are  5  to  7  days.  Along  the  mainland  shore  of  the  gulf, 
the  waves  lose  some  energy  due  to  topographic  irregularities  of  the  main  basins  of  the  gulf  After  the 
waves  reach  the  sill  (at  y  ~  900  km  in  Figure  2)  where  the  depth  decreases  abruptly  from  1500m  to 
200m,  a  small  fraction  of  the  energy  continues  north  where  it  is  dissipated  in  the  shallower  water.  Most 
of  the  energy  is  reflected  at  the  sill  as  a  wave  that  travels  south  along  the  coast  of  the  Baja  Peninsula. 
This  process  is  illustrated  in  Figure  2.  The  wave  that  leaves  the  gulf  on  the  Baja  side  contains  only  a 
small  fraction  of  the  incident  energy. 

Wind  stress:  The  dominant  winds  are  aligned  along  the  gulf  and  are  southward.  Wind  reversals  occur, 
however,  so  that  both  sides  of  the  gulf  experience  some  upwelling  favorable  conditions.  The  yearly- 
averaged  wind-forced  surface  currents  flow  out  of  the  gulf,  while  the  average  flow  at  depth  is  into  the 
gulf  This  feature  is  consistent  with  observations  (Bray,  1988)  and  helps  explain  the  high  productivity  in 
the  gulf  The  most  striking  feature  of  the  wind-forced  circulation  is  the  formation  of  eddies  in  the  south 
and  central  gulf  The  scale  of  these  eddies  is  about  90  km  and  they  are  generated  on  both  sides  of  the 
gulf  Numerical  experiments  without  surface  heat  flux  show  that  the  intensity  of  the  eddies  is 
substantially  increased  by  the  presence  of  the  strong  surface  heat  flux  in  the  gulf  (annual  average  90 
W/m^). 

Tides:  The  external  M2  tide  is  sufficiently  energetic  (2m  and  1.5  m  s  ')  to  generate  a  strong  residual 
circulation  (max  velocities  0.18  m  s  ')  in  the  vicinity  of  the  sill  and  in  the  region  farther  north.  There  is 
substantial  generation  of  internal  tides  (20-30  m  amplitude)  in  the  neighborhood  of  the  sill.  Mixing 
produced  by  tidal  forcing  is  important  in  the  sill  region.  The  total  energy  conversion  as  a  result  of  tidal 
mixing  is  comparable  to  the  total  energy  dissipated  by  bottom  friction. 


IMPACT/APPLICATIONS 


The  numerical  modeling  studies  with  the  primitive  equations  applied  to  both  wind-forced  inner  shelf 
flows  and  to  wave-averaged  currents  in  the  surf  zone  constitute  important  initial  steps  toward 
development  of  a  single,  unified  modeling  approach  to  three-dimensional  circulation  processes  in  the 
nearshore  inner  shelf  region.  Modeling  studies  have  provided  new  quantitative  information  on  the 
nature  of  the  wind-  and  coastal- trapped  wave-driven  mesoscale  circulation  in  the  Gulf  of  California. 

TRANSITIONS 

RELATED  PROJECTS 

Some  aspects  of  the  primitive  equation  Princeton  Ocean  Model  studies  of  surf  zone  flow  fields  are 
jointly  funded  by  ONR  Grant  N00014-99-1-1051,  (NOPP)  "Development  and  Verification  of  a 
Comprehensive  Community  Model  for  Physical  Processes  in  the  Nearshore  Ocean". 
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Figure  1.  Results  from  two-dimensional  POM  simulations  of  shelf flows  at  Duck,  NC.  The  two 
upper  panels  show  depth-averaged  alongshore  currents  at  13  m  and  26  m  from  both  the  data  (blue) 
and  the  model  (red)  for  the  period  August  7, 1994  to  November  11, 1994.  The  corresponding 
correlation  coefficients  are  also  included.  The  six  lower  panels  give  monthly  mean  fields  of  the 
potential  density  (Oe),  streamfunction,  and  alongshore  velocity  and  their  standard  deviations  for 
August  (stratified  shelf)  and  October  (unstratified  shelf).  Standard  deviations  are  in  color,  while 

the  monthly  means  are  overlaid  as  black  lines. 
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Figure  2.  Results  from  numerical  model  studies  of  the  mesoscale  circulation  in  the  Gulf  of 
California.  This  experiment  involves  forcing  by  coastal-trapped  waves  that  propagate  northward 
into  the  gulf  along  the  coast  of  mainland  Mexico.  The  propagation  and  decay  in  the  gulf  of  a 
coastal-trapped  wave  disturbance  are  evident  in  the  sea  level  elevation  (color  contours)  and  depth- 

averaged  velocity  vector  fields  from  8-14  July  1984. 


